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Penetration and energy-loss theory of electrons in solid targets 

K KANAYAt and S OKA YAM . 

Elect rotechnical Laboratory, Mukodai-c no 5-4-1 , Tanashi-shi, Tokyo, Japan 
MS received 5 April 1971 , in final form 4 October 19?! 

Abstract. Starting from a simple atomic model giving the potential between electrons 
and atoms as V{r)-Ze 2 a 9 - X jsr t with the empirical value j=fl, we combine the diffusion 
effect due to multiple collisions and the energy retardation in accordance with a modified 
Thomson-Whiddington law, with the scattering cross section in the Lenard absorption 
law. On this basis, consistent expressions are obtained for the fraction of transmitted 
electrons in amorphous solid targets, the backscattering fraction with depth, the fraction 
of electrons absorbed per unit mass-thickness and the depth-dose function, which arc 
in good agreement with experiments over the energy range 10- J 000 kcV. 

A diffusion model represented by a sphere whose centre is located at the maximum 
energy dissipation depth, related to the diffusion depth and the range, is found to 
agree well with experiments, 

1. Introduction 

An explanation of the phenomena connected with electron penetration into solid 
materials requires information about the attenuation of the electron beam in both numbers 
and energy. Electron penetration into solid materials has been studied for over 70 years. 
Most of the relationships have been obtained either by experiments with a high-energy 
beam (Lenard 1 895) or by theoretical considerations under assumptions which are also 
applicable only at high energies (Bothe 1933, 1949, Bethe 1933, Bethe et al 1938). 
Furthermore, important papers by Wentzel (1922, 1927), Moliere (1947, 1948), Lenz 
(1954), Marton et al (1962) and Smith and Burge (1963) were concerned with plural 
scattering events. The results have been presented in a collected and generalized form by 
Cosslett and Thomas (1964a,b, 1965), based on the experimental data with the quantita- 
tive discussions. Many experimental results have been reported in the field and compared 
with the different collision models and various theories by Goudsmit and Saunderson 
(1940), Spencer (1955), Lenz (1958), Meister (1958), Archard (1961), Everhart (1960), 
Makhov (1960a,b,c), Shimizu and Shinoda (1963) and Dashen (1964). 

A simple and semi-empirical theory of electron scattering in solid materials may be 
required in connection with electron probe microanalysis, scanning electron microscopy, 
electron beam microrecording and micromachining. When a stream of electrons 
penetrates into a solid target, electrons maybe scattered either elastically or inelastically. 
Electronic stopping is due to inelastic collisions with atomic electrons in which the 
incident electron excites or ejects atomic electrons with loss of energy. The corresponding 
momentum transfer is small because electrons are light particles, but the energy loss is 
very large. Nuclear stopping, on the other hand, arises from nearly elastic collisions 
with atomic nuclei, with transfer of both energy and momentum. 

t Now at University of Kogakuin, Nishi-shinjuku 1-24-2, Shinjuku-ku, Tokyo, Japan. 
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Thus electrons are supposed to travel straight into the target, suffering energy loss 
due to the electronic collisions, and also to be deflected by the nuclear collisions. Back- 
scattering electrons are those elastically scattered mainly due to nuclear scattering. 
The plasmon excitation also affects the energy loss of electrons in the solid. Even though 
the cross section of plasmon excitation is greater, the energy loss seems to be generally 
small compared with the energy dissipation due to the electronic collision. Therefore 
the plasmon effects are neglected in this simple theory. 

A quasi-elastic scattering theory by Lindhard et ai (1963) for a stream of ions 
penetrating into a solid target, which takes into account both types of stopping, may be 
applied. The cross sections of the nuclear and electronic collisions are derived from an 
atomic model represented by the power potential V(r) = Ze 2 a 8 ~ 1 /sr s y with j= J, and 
the maximum range for the accelerating voltage from 10 to 1000 keV, which is pro- 
portional to the incident energy Eo w , is compared with the most reliable calculation by 
Bergcr and Seltzer (1964) and the experimental results Cosslett and Thomas (1964b). 

Accordingly, the semi-empirical expressions for the fractions of transmission, back- 
scattering and absorbed energy are obtained as a function of the reduced depth y — x/R, 
which may be available for different target materials and the wide energy range. 
Furthermore, after comparing the results obtained by the present theory with the Archard 
diffusion model, we propose a modified diffusion model. The model, represented by a 
sphere whose centre is located at the maximum energy dissipation depth, related to the 
diffusion depth and the range, may give an explanation to the experiments by many- 
authors. 

2. Nuclear and electronic scattering cross section 

After Lindhard's theory concerning ion-beam scattering, the potential K(r) between the 
electron and the target atom is assumed to be 

V(r)=Ze 2 a'- l /sr* (I) 

with 

a = 0-8853a H Z-i/3 

where Z is the atomic number of the target, e the electronic charge and a the effective 
screened radius of the atom. The number 0-8853 = (97r 2 )i/32 -7/3 j s t fc e familiar Thomas- 
Fermi constant and an is the Bohr radius of the hydrogen atom (5-292 x 10~ 9 cm), s is 
a numerical parameter: .y=l corresponds to Rutherford scattering and ^ = 2 to 
d£/dx = constant. The differential scattering cross section for angular scattering is 
given by 

d ^ = 273/ 2 /«< /1/2 ) d ' (2) 

with 

f a ( t m) = x tf 

where r= e 2 sin 2 <f>/2; t = a/b and <f> is the deflection angle in the centre of gravity system. 
b is the collision diameter and A, is a constant determined empirically. On the other 
hand, the differential scattering cross section for the energy transfer is derived from 
(2) and t = {a/b)*T/T m : 

770 2 /b\ 2 " T m u> ... 
dar=A._y j^dT (3) 

where T and T m are the energy transfer and its maximum value. 
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For the special case of J= J, which gives good agreement with empirical results 
in the energy range, the differential nuclear scattering cross section is rewritten from (2) 
and& = Ze 2 /E: 

. 2"/«7rai/3<>io/3 Z 5/s s in0d0 u 

d ^- A * 4 £5/3 (p+cos 6 yv* { ' 

where 

6=n — <f>. 

When we consider the electron penetration in the target, electrons suffering deflections 
more than 90° do not travel into the subsequent layers of the target and the cross section 
is given by integration from 0 = 0° to 90° as a first approximation. As was pointed out 
by Archard (1961), many electrons are deflected between 0° and 90°, and some of them 
arc lost by a multiple collision effect which we will now describe. Consider en electron 
initially deflected at 45 ° : it must suffer a second deflection of 45 ° in the opposite sense 
to get back into its original direction ; but it might equally be deflected by 45 ° in the 
same sense, thereby acquiring a total deflection of 90° and becoming lost 10 subsequent 
layers. Thus a second approximation for the cross section has been given by adding 
half the integration from #--=90° to 135°. But the correction may be not enough, since 
the triple and more collisions must exist in fact, as shown in figure 1 . The better approxi- 
mation may be expressed by 

J 0 (l +CO~S '0)11/6 +i J n/2 (1+COS0)""/ 6 (1+ COS 0)11/0 

and the value of 3-33 is used as an empirical value in the following calculation. 



Incident beam 




Figure 1. Modified diffusion model of electron-beam penetration in a target: R is the 
maximum range; xd the diffusion depth; xg the maximum energy dissipation depth; 
tb the backscattering range; tan flo-na/xE. (1), (2) and (3) refer to the number of 
times the electrons are deflected. 
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Then the total scattering cross section for the angular deflection due to the nuclear 
collisions is expressed by 

, 2-977rfll/3 e 10/3 Z S/3 ™ 

* n=A ' iw • (5) 

On the other hand, the differential cross section for energy loss due to the ctronic 
collisions is derived from (3) and 6 = 2e 2 /£: 

d. e =A s ^. 2 / 2 ^ 5/3r ™ 5 . /6dr (6) 
e ' 2 \aE) r»/« 

and resultant energy loss is obtained as follows : 



dx Jo dT 



__ , 3 X 2 5 / 3 7rfll/3 e 10/37VZ n , 
~ Aa £2/3 " K) 

where N is the number of atoms per unit volume in the target. The two equations 
(5) and (7) are the fundamental relationships of clastic (nuclear) and inelastic (electronic) 
scattering theory. 



3. Range-energy relationship 

The maximum range can be derived from the energy-loss equation (7) : 
dE 



ZT-5/3 

(8) 



dE/dx 
£o 5/3 



X s x 5 x 25/3 7ra i/3 e io/3 jVZ 

where the number N is given by N=N & p/A ; N & is the Avogadro number, 6-023 x 10 23 ; 
P (g cm -3 ) the density; A (g) the atomic weight; and Eq (eV) the incident energy. Then 
the mass-range pR can be expressed by 

P* = 5-025 x 10" 12 /l£o 5/3 /A a ZS/9. (9) 

This relationship satisfies qualitatively the experimental results of Makhov (i960a,b,c) 
and Schumacher (1967) in respect of the dependence of incident energy Eo and atomic 
number Z. Good agreement with experiments from 10 to 1000 keV is obtained by 
taking A A = 0*182. Equation (9) can be rewritten under the relativistic corrections of 
energy £ 0 (1 +0 978 x 10-«£ 0 )/(l + 1-957 x W' 6 E 0 ) and a=0-8853* H Z-i/ 3 x (1 + 1-957 x 
lO- 6 ^)- 1 : 

2-76xlQ-iME 0 5/3 (i +Q-978 xW -* Ep)W 
P Z**» (1 + 1-951 x\0-*Eq)W ' 

Figure 2 shows a comparison of the voltage dependence of mass-range for several 
target materials by (10) with the experimental (Cosslett and Thomas 1964b) and calcu- 
lation results (Berger and Seltzer 1964). 
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From equations (7) and (8) the energy E of electrons at depth x is simply expressed 
in terms of the reduced depth y — x/R: 

Using the above relation of energy retardation, we can obtain, consistently, the fractions 
of transmission, backscattering and absorption as functions of y. 



4. Transmission and diffusion depth 

The general form of variation with thickness in the fraction t)t of the incident current 
which is transmitted into the forward hemisphere as shown in figure 1 may obey an 
exponential relation similar to the Lenard law — 

i?T = -=exp(- N<TnX). (1 2) 

The total scattering cross section a n (equation (5)) is connected with the fractional 
range-energy relationship (11) such that 

, 2-977Tfll/3 5 10/3 Z 5/3 
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By substituting (8) and (13) into (12), t> t can be expressed as a function of y and y: 

rn-exp^-^J. (14) 

The parameter of y involves the effects of diffusion loss due to multiple collisions for 
return^ electrons and energy retardation due to electronic collisions. It is related to 
the atomic number alone as follows: 

y — ^3 Z 2/3 = 0*I87Z 2/3. ( 1 5) 

5 x 2 5 ' J 




Figure 3 shows the variation of r} T in several target materials as a function of y. For an 
clement of low Z (such as carbon), the value of tj T is too small in the region 0^j><0-3 
because of the overestimation of scattering cross section in this region. For elements of 
high Z (such as copper, silver, gold and uranium), the values calculated by (14) agree 
closely with the experimental results by Seliger (1955) and Cosslett and Thomas (1964a). 
Fractional electron transmission 7jx for copper as a function of mass-thickness px is 
shown in figure 4. 

From the definition of the diffusion depth xry where the transmitted fraction is 1/e, 
it follows that 

**- f — — (16) 

which seems to be more reasonable than the expressions given by Meister ( 1 958), Archard 
(1961) and TomJin (1963) and accords well with the Monte Carlo calculations by Bishop 
(1965, 1967). The calculation results from (16) are compared with the experiments by 
Cosslett (1964) and Seliger (1955) in figure 5. 
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px f/jq cur 1 ) 

Figure 4. Fractional electron transmission tjt for Cu (29) (y = 1 *76) as a function of mass- 
thickness at several values of incident energy £o, compared with experimental results (#) 
of Cosslctt and Thomas (1964a). 




5. Backscattering 

The general form of the variation with thickness in the fraction ?;b of the backscattered 
electrons, which are deflected inside the limiting angle 6 0 subtended by the backscattered 
radius rn of the centre of a sphere model, as in figure 1, is assumed to have the same 
exponential form as tj T , but the absorption factor yb must be larger than y because of 
diffusion loss due to multiple collisions. 
5 
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By differentiating (14) with respect to >«, using (4) and normalizing by 

1 



n = 



J* sin 0d 0/(1+ cos 3) n ' 6 



we obtain 



dy" ' 



' dy 



f a ° si 

« I 

Jo 



sin 0d0 



cos 6) n '* 



where £>(» is the diffusion-loss function of backscattered electrons with respect to y. 
It may be simply expressed by , 



Ay 



sin 0 d0 



<l-,)« CXp i l-WJo (l+< 



cos 0) 11 ' 6 
sin 0d0 



cos 0) 11 ' 6 



with 



YB = 1 *9y. 

Then the backscattering fraction rj B can be obtained by integration: 

sin 6 d0 



(17) 
(18) 



With cos 0o= y /(I — ,y), this leads to 



cos S)W 



The fractional backscattering 773 from several target materials as a function of reduced 
depth y, as calculated by (19), is shown in figure 6. The present result is rather in closer 




Figure 6. Fractional backseat terings rj B from several targets as a function of reJuced 
depth y y compared with experimental results of Cosslett and Thomas (1965) for Cu (©) 
and Au (O) at 5-20 keV. 
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agreement with experiment (Cosslett and Thomas 1965) than the existing theories 
(Archard 1960, Everhart I960, Archard and Mulvey 1963). 



6. Absorption of electrons 

Following Cosslett and Thomas, we obtain the fraction rj A of the incident beam absorbed 
in a solid target by difference: 

VA = 1 — (tjt+^b). 

If rj T is substituted by (14) and t?b by (19) into the above, rj A becomes 

The absorbed fraction calculated by (20) is in good accord with the results of Cosslett 
and Thomas (1965), as shown in figure 7. 




Figure 7. Relative proportions of electrons transmitted (i?t), backscatterci (^b) and 
absorbed (*?a) in Cu as a function of the reduced depth y. Broken curves and points 
indicate experimental results of Cosslett and Thomas (1965) for Cu at JOkeV. 



The distributions of absorbed electrons at different incident energies can be repre- 
sented by a single curve for a given element when the normalized fraction is plotted 
against depth as a fraction of the range y 9 as suggested by Cosslett and Thomas (196 j). 
The normalized fraction of electrons absorbed per unit mass-thickness is obtained by 
differentiating (20): 

x(^ 6 -(l->0 5 >*)}. (21) 
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x!0 3 6| 1 1 \ 1 1 1 1— 1 1 




Figure 8. Fraction of electrons absorbed per unit mass-thickness as a function of depth 
px in solid targets of Cu and Au at 5, 10, 15 and 20 keV. The values at 5 keV arc drawn 
to half scale. 



The fraction of electrons absorbed per unit mass-thickness as a function of depth .v 
in solid targets of copper and gold at 5, 10, 15 and 20 keV is shown in figure 8. The 
distributions for the two elements are very similar at the same incident energy. The 
variation shows qualitative agreement with the differential absorption evaluated from 
experimental data of transmitted and backscattered electrons by Cosslett and Thomas. 

7. Energy dissipation 

The dissipation of energy with depth may be calculated in a similar way to the treatment 
by Cosslett and Thomas (1965): 

£a = Eq — VT E — rjB Eb 

where E and Eb refer to the energy of transmitted electrons through matter, reduced by 
electronic collisions as given by (1 1), and the mean energy of backscattered electrons from 
a target due to nuclear collisions. Values of Eb/Eo are provided by several authors 
(Brand 1936, KulenkampfT and Riittiger 1954, KulenkampfT and Spyra 1954, Sternglass 
1954, Klein 1968), as shown in figure 9. The value that best fitted the curve shown 
in the figure was used for the present calculation. 

We deduce that, corresponding to the above equation, the energy E\ absorbed in 
the fractional layer of material between the surface and depth x is given by 
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Mean fractional energy of reflected electrons En/£o against Z, the atomic 
Experimental points: O (Brand 1936 — 32 keV, 5 = 8°); • (Sternglass 1954 — 
); □ (Kanter 1957— lOkcV, 0=15°); B (Kantcr 1957— 10-70 keV, 




Figure 10. Fraction of energy dissipated £,\/£o within a given fraction of the range v 
in targets of C, Al, Cu, Ag, Au and U. 



Figure 10 shows the result calculated by (22) for several target materials. The calculated 
distributions of the energy explain fairly well the fact that at a given incident energy the 
curves are first almost coincident, for the energy is retarded due to electronic collisions, but 
then diverge increasingly owing to the difference in the amount of energy lost by back- 
scattering due to nuclear collisions. At the end of the range this is equal to the 
amount [i7b]*-i/2 (£b/£o). 
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The reduced fraction of energy dissipated in unit mass-thickness d(£A/Eo)/d(/>x) 
can be obtained by differentiating (22): 

(23) 

2 5 — — 



20 
1-5 

— H 
^ -o 10 
cx 

0 5 



u 0 1 0-2 0-3 0-4 0-5 0 6 0 7 0 8 0 9 10 

/ 

Figure 11. Normalized fraction of energy dissipated in unit mass-thickness as a function 
of the reduced depth y. 



Figure 1 1 shows the result calculated by (23) for several target materials. In the figure 
the peak position corresponds to the fraction of maximum energy dissipation depth. 
Its values are 0-56, 0-38, 0-24, 0-18, 014 and 0-12 for carbon, aluminium, copper, silver, 
gold and uranium. Figure 12 shows the depth-dose function per unit mass-thickness, 
d£ A /d(px), for two typical target materials of copper and gold, compared with the results 
by Gosslett and Thomas (1965). The theoretical and experimental distributions show 
only qualitative agreement, but the peak value is located at almost the same reduced 
depth. Better agreement would be possible by selecting more appropriate values of 
pR and E^jEo, 

8. Diffusion model 

On the basis of the Bethe energy-loss theory, Archard (1961) proposed the so-called 
'diffusion model' of a sphere in which electrons move equally in all directions from the 
depth of complete diffusion xj> in such a way that their overall paths are equal to the 
full range R. For a high atomic number this model agrees fairly well with the photograph 
of electron glow published by Ehrenberg et at (1953, 1963) ; but for a low atomic number 
it cannot be applied, because backscattered electrons reaching the surface in their 
original directions vanish. This is because the model disregards electrons undergoing 
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Figure 12. Depth-dose function in solid targets of Cu and Au at 5, 10, 15 and 20 kcV, 
compared with experimental results of Cosslctt and Thomas (1965) for Cu (•) at 
15 kcV and for Au (O) at 10 keV. 



large-angle elastic reflection between the surface and the depth of complete diffusion. 
Accordingly, a modified diffusion model may be assumed which represents the distri- 
bution of electrons by a sphere whose centre is located at the depth xe of energy dissi- 
pation; this takes account of both the absorbed electrons and deflected electrons. 
The backscattering range rs is given by 

r B =^ (24) 
l+y 

where the best fitting is obtained by taking C=l-I. 

The maximum energy dissipation depth x E can be obtained by a simple geometrical 
relation : 

* E= 2 l^oT^J 2(i+ y )V • (25) 

The values of xe/R agree with those for the peak in the fractional energy dissipation 
curves (figure 11). The values of xe/R are 0-41 (for carbon), 0-34 (aluminium), 0-25 
(copper), 0-20 (silver), CM 4 (gold) and 01 2 (uranium), compared, with the positions of 
the energy dissipation peaks : 0-56, 0/38, 0-24, 0 18, 0-14 and 012. In elements of low Z, 
even though the difference is greater than for elements of high Z, we can assume as a 
first approximation that electrons spread out from points which coincide with the 
maximum energy dissipation depths. 

Then the backscattering coefficient r can be obtained by the following equation: 



1 fV 



2rr sin 6 d0 = K! - cos 0 O ) (26) 

H7T J q 

with 
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The backscattering coefficient r from (26) accords well with the value nu(>— i) found 
from (19), as shown in figure 13. The values of r are 015 (for carbon), 0-24 (aluminium), 
0 33 (copper), 0-38 (silver), 0-42 (gold) and 0-43 (uranium), compared with the back- 
scattering fractions -q B at the end of the range 0> = i): 018, 0-25, 0-33, 0-37, 0 ; 40and 0-41. 
Figure 14 represents the modified diffusion models for several target materials, where the 
broken curves][are the original diffusion models. 
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Figure 13. Backseat tering coefficient r for atomic number Z according to tlr* modified 
diffusion model, compared with the results of different authors. Broken curve indicates 
the backscatiering fraction 77s at the end of the range 0=i). Experimental points were 
collected by Archard (1961); references will be found therein. 

9. Conclusions 

The results of interest in connection with the fundamental theory of electron scattering, 
such as mass-range, transmission, diffusion depth, backscattering, energy loss and 
maximum energy-loss depth, are consistently expressed in normalized form as a function 
of the reduced depth y. And, based on the relationship existing between backscattering 
and energy dissipation, a diffusion model represented by a sphere whose centre is located 
at the maximum energy dissipation depth is proposed. 

The present treatment gives only a first-order approximation to the scattering events 
of penetrating electrons in a solid target. Better quantitative comparison between 
theoretical and experimental results may be obtained by using corrected values of the 
parameters A,, y and which have been arbitrarily assumed on an empirical basis 
because they are independent of the fundamental course of the theory. However, in 
the case of high voltages above 1 MeV, the power-potential law with j = f cannot be 
applied. It seems that the Rutherford scattering theory and Betbc energy-loss law are 
more realistic expressions in such a case. According to the same treatment, by using the 
Rutherford scattering model with the absorption parameter y is given by 
cZ 

y ~lnOEo/I3-5Z) 

where 

c =0 1-0-3. 
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Figure 14. Representation of modified diffusion models for several targets: 
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The distributions of tjt, >?b, t?a and Ea vary slightly with the accelerating voltage. If 
more reliable data for comparison become available, the continued treatment by the 
Rutherford scattering model will be reported elsewhere. 
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